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METODOLOGÍAS PARA LA SIMULACIÓN 
DE LA FLUIDODINÁMICA EN UN MOTOR 
DE COMBUSTIÓN INTERNA MEDIANTE 
SOLIDWORKS 


Methodologys for the simulation of fluid 
dynamics in an internal combustion engine 
through SolidWorks 


RESUMEN 


La representación computacional de los efectos fluidodinámicos durante la interacción entre el aire como 
fluido de trabajo y las diferentes geometrías que conforman las carreras de admisión y escape en un Motor 
de Combustión Interna (MCI) mediante simulaciones de Dinámica de Fluidos Computacional (CFD) 
usando inéditamente SolidWorks Flow Simulation (SWFS) es reportada en este trabajo. Para lo cual se 
necesitó comprender los requerimientos del proceso de configuración de SWFS, adaptar las condiciones 
iniciales y de frontera propuestas en antecedentes bibliográficos y usar una geometría del MCI albergada 
en una base de datos geométrica de acceso libre. En los resultados fue posible exponer las características 
del fluido de trabajo al interior del cilindro para varias aperturas de válvula de admisión y escape en términos 
de variables claves como presión, flujo másico y campos de velocidad denominados Swirl y Tumble; 
mediante la validación fue posible evidenciar que las metodologías propuestas representan la 
fluidodinámica del aire con un alto nivel de convergencia según trabajos reportados en la literatura. 


Palabras clave: Dinámica de fluidos computacional; Ecuaciones de Navier Stokes; Motores de combustión 
interna alternativos; Carrera de admisión; Carrera de escape. 


ABSTRACT 


The computational representation of the fluid dynamic effects during the interaction between air as a 
working fluid and the different geometries that make up the intake and exhaust strokes in an Internal 
Combustion Engine (ICE) by means of Computational Fluid Dynamics (CFD) simulations using previously 
unpublished SolidWorks Flow Simulation (SWFS) is reported in this work. For which it was necessary to 
understand the requirements of the SWFS configuration process, adapt the initial and boundary conditions 
proposed in bibliographic precedents, and use an ICE geometry housed in a freely accessible geometric 
database. In the results it was possible to expose the characteristics of the working fluid inside the cylinder 
for various intake and exhaust valve openings in terms of key variables such as pressure, mass flow and 
speed fields called Swirl and Tumble; through validation it was possible to show that the proposed 
methodologies represent the fluid dynamics of air with a high level of convergence according to works 
reported in the literature. 


Keywords: computational fluid dynamic; Navier Stokes equations; alternative internal combustion engine; 
intake stroke; exhaust stroke. 
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INTRODUCCIÓN 


Los motores de combustión interna (MCI) cumplen un papel crucial en diversas aplicaciones tales como 
transporte y generación de energía para la sociedad. Sin embargo, estas máquinas están ligadas al uso 
tradicional de combustibles fósiles (FF) generadores de gases nocivos tales como los óxidos de nitrógeno 
(NO,), monóxido de carbono (CO) entre otros gases contaminantes, que en conjunto con las ineficiencias 
asociadas a las pérdidas térmicas, mecánicas y de flujo que presentan estos equipos en su funcionamiento 
normal; reflejan la necesidad de actividades de investigación y desarrollo (I+D) que busquen un modo de 
operación eficiente y sostenible de los MCI (Lalwani et al., 2019). 


Técnicamente, un MCI es un conjunto de componentes mecánicos que permiten transformar la energía 
química de un combustible líquido o gaseoso en energía térmica mediante la combustión y a través de la 
expansión del fluido de trabajo (mezcla aire — combustible), convierte la energía térmica en trabajo mecánico 
(Jacobs, 2020). En el funcionamiento normal de los MCI existen algunos factores previos a la combustión, 
los cuales influyen directamente en las prestaciones, la eficiencia y la formación de emisiones contaminantes 
del motor, dichos factores se dan durante el proceso de formación del fluido de trabajo y son los patrones 
de rociado del combustible, el movimiento del fluido de trabajo en el cilindro y las características físicas de 
la mezcla (Chinnamuthu et al., 2021); debido a lo anterior, investigadores aseguran que una forma de 
abordar el estudio y comprender la fenomenología de estos factores claves es mediante la dinámica de 
fluidos computacional (CFD) (Hsieh et al., 2021). Esta técnica de estudio combina la mecánica de fluidos 
moderna, con métodos numéricos e informática, por lo que una simulación CFD bajo ciertas condiciones y 
parámetros dados es equivalente a un experimento numérico, escatimando de esta forma tiempo y recursos 
ante alternativas analíticas y experimentales (Liu, 2021). 


En la literatura hay diversos trabajos que reportan el uso de la CFD en el análisis de las características 
fluidodinámicas previas a la combustión en los MCI mediante muchos programas que integran esta técnica 
de estudio y según la revisión bibliográfica, SolidWorks Flow Simulation (SWFS) ha sido poco explorado en 
este campo. Esta herramienta CFD es un complemento integrado al software de diseño asistido por 
computadora SolidWorks, la cual permite obtener soluciones parciales de las leyes de conservación de masa, 
momento y energía, conocidas también como las ecuaciones de Navier Stokes, que en conjunto con la 
ecuación de continuidad gobiernan el movimiento de los fluidos (Matsson, 2021). Los pocos trabajos 
existentes con SWFS se enfocan en el estudio de cómo influye el diseño del colector de admisión en las 
características del fluido de trabajo y las prestaciones de un MCI diesel transformado a gas licuado de 
petróleo (GLP) (Jemni et al., 2011); el enriquecimiento de Hidrógeno (H2) y su impacto en las 
características fluidodinámicas en el cilindro, las prestaciones del MCI y las emisiones de escape de un motor 
que normalmente funciona con GLP (Kacem et al., 2016); el análisis de las características del fluido de 
trabajo en el cilindro con estrategias innovadoras de admisión para un MCI multiválvulas (Nigro et al., 2017) 
y en la influencia de la ubicación de la inyección de H, en las características fluidodinámicas, el rendimiento 
y las emisiones de un MCl que funciona con GLP (Jemni et al., 2018). 


Cabe resaltar que la CFD no solo aborda el estudio de los MCI, se expande a muchos escenarios y entre 
ellos la medicina cardiovascular, donde se usa como método /n silico no invasivo para caracterizar flujos 
sanguíneos en el sistema cardiovascular en condiciones fisiológicas o patológicas con aplicaciones 
potenciales en la toma de decisiones clínicas, planificación, estrategias de terapia, diseño y optimización de 
dispositivos médicos (Totorean et al., 2022). 
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En este orden de ideas, al considerar que la técnica CFD es una herramienta poderosa que puede ser usada 
en muchos escenarios de interés para las distintas comunidades científicas y que el uso de SWFS, según la 
revisión bibliográfica, es poco reportado en estudios sobre la interacción entre el fluido de trabajo y las 
diversas geometrías que conforman a un MCI. El objetivo de este trabajo apuntó a la representación 
computacional de la fluidodinámica presente en la interacción del aire como fluido de trabajo durante las 
carreras de admisión y escape con las geometrías que conforman un MCl mediante una simulación CFD en 
el entorno de SWFS y la documentación del proceso tras dicha representación, el cual incluyó la identificación 
de las condiciones iniciales y de frontera, la adaptación de estas a los requerimientos de SWFS, en donde 
destacó un diferencial de presión para hacer fluir el aire a través de las geometrías del motor; el análisis de 
mallado y el postprocesamiento de resultados en función de variables como flujo másico, presión y los 
campos de velocidad (Swirl y Tumble); la validación de resultados para los procesos de admisión y escape 
se realizó con trabajos afines reportados en la literatura. 


METODOLOGÍA 


Para la consecución del objetivo planteado en este trabajo, fue necesario comprender los requerimientos 
tras el proceso de configuración de una simulación CFD en el entorno de SWFS, el cual se compone por las 
etapas básicas de preproceso, procesamiento y postproceso, ilustradas en la Figura 1. 


PROCESAMIENTO 
PREPROCESO CS men dcir aca 
= Verficar 
1. Modelo CAD 1 Create Lips/Saliente extruir geometría | 
1. Sistema de unidades. 1 +— + ' 
2. Ti terísticas fisi . 
del ca dr si . 1. Dominio computacional I POSTPROCESO 
2.1. Análisis interno o externo 1 2. Subdominio del fluido Il 4.Malla 
2.2. Conducción de calor, radiación, 1 3. Condiciones de frontera l 2. Gráfico de contornos 
dependencia del tiempo, gravedad, y 1 [Entradas/Salidas] ' 3. Gráfico de superficie 
rotación. 3.1. Condiciones de flujo ' 


3. Tipo y características del 
fluido. 


3.1. Gases, líquidos, no newtonianos, 


compresibles, gases reales y vapor. 
3.2. Régimen laminar, turbulento y 
ambos. 
3.3. Cavitación. 
4. Condiciones de pared. 


4.1. Pared adiabática, flux de calor, 


tasa de transferencia de calor y 
temperatura en pared. 
4.2. Rugosidad. 
5. Condiciones iniciales 
[Parámetros] 
5.1. Termodinámicos [T;, P;]. 
5.2. Velocidad [X, y, Z]. 


5.3. Turbulencia [Intensidad y Duración 
(1 — L)/Energía y Disipación (K — e)] 


5.4. Concentración [fracción de 
volumen o masa] 


- Velocidad, flujo másico y 
volumétrico 
3.2. Condiciones de presión 
- Ambiente, estática y total 
3.3. Condiciones de paredes 
- Ideales y reales 
4. Objetivos de simulación 
[Parámetros] 
4.1. Generales 
4.2. Punto especifico 
4.3. De superficies 
4.4. Volumen 
4.5. Insertar ecuación de objetivo 
5. Mallado 
5.1. Malla general 
5.2. Malla local 
5.3. Crear malla 


InputData 


CFD en SWFS. 


4. Isosuperficie 

5. Trayectoria de flujos 

6. Estudio de partículas 
6.1. Configuración 

7. Parámetros de punto 
especifico 

8. Parámetros de superficie 
8.1. Exportar a Excel 

9. Parámetros de volumen 
10. Gráfico XY 

11. Gráfico de objetivos de 
simulación 

12. Reportes 

13. Animación 

14. Exportar resultados 


Figura 1. Requerimientos generales de las etapas básicas de preproceso, procesamiento y postproceso para una simulación 
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Este diagrama de flujo general se estableció debido a la revisión de trabajos como el de Rojas $ Huaraz 
(2018) y Tordini (2019), manuales técnicos de SWFS (Dassault Systemes, 2021) y la interacción que se 
tuvo con el software SWFS, lo cual permitió identificar algunos requerimientos claves al momento de 
configurar una simulación CFD tales como la selección del modelo o diseño geométrico asistido por 
computadora (CAD), el establecimiento de las condiciones iniciales y de frontera, tipo de análisis, tipo y 
características del fluido de trabajo, configuración de mallado y la selección de los objetivos de la simulación; 
los cuales se tratan de forma específica en las siguientes subsecciones. 


Modelo CAD 


El modelo CAD seleccionado para este estudio se ¡ilustra en la Figura 2, este cuenta con los siguientes 
componentes de un MCI: múltiples de admisión/escape, válvulas de admisión/escape, cilindro y pistón. Este 
modelo se encuentra disponible para fines académicos en la base de datos (GrabCAD, 2021). 


-s- 0 


Figura 2. Modelo CAD ilustrado en las vistas (A) derecha, (B) isométrica, (C) superior, (D) frontal (admisión) y (E) posterior 
E 


Condiciones iniciales y de frontera 


Para estos requerimientos se tomó como punto de partida las condiciones reportadas en el trabajo de 
(Scaggion, 2018) encontrado en la revisión bibliográfica, las condiciones iniciales y los objetivos para el 
proceso de admisión se listan en la Tabla 1. 
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Tabla 1. Condiciones iniciales (variables independientes) y objetivos (variables dependientes) del proceso de admisión donde 
AP es un diferencial de presión, haqm apertura de válvula de admisión y T la temperatura para todo el sistema. 


Proceso de admisión 


Condiciones iniciales Objetivos de simulación 
) 
AP = 5.000 Pa á ; 
(con respecto a la presión atmosférica) Campos de velocidad del fluido 
hadm =1,2,3,4,5,6 y 7 mm Campos de presión 
T =313K Flujo másico 


Debido a que interesa conocer como es el comportamiento fluidodinámico del fluido de trabajo al interior 
del motor, específicamente en el cilindro, se tienen en cuenta las siguientes consideraciones para las 
simulaciones: el fluido de trabajo es aire, no se intercambia calor con el entorno y las paredes se consideran 
con rugosidad cero, la presión atmosférica se toma como 1 atm; en cuanto a las condiciones iniciales de 
turbulencia y concentración, se trabaja con la configuración predeterminada de SWFS. 


Para la simulación del proceso de escape, se tomó en cuenta como condición de frontera en el múltiple de 
admisión la magnitud del flujo másico admitido cuando hagm= 4 mm y la presión atmosférica a la salida del 
colector de escape; la válvula de admisión tuvo una apertura constante de 4 mm mientras que la de escape 
(hesc) variará según el rango de aperturas hese = 1,4 y 10 mm; las variables dependientes para este 
escenario serán flujo másico a la salida del colector de escape, velocidad de admisión y de escape, campos 
de velocidad y presión. 


Análisis de mallado 


SWFS integra la opción de malla automática, la cual permite definir su escala entre 1, que representa el 
nivel de malla con celdas más grandes y 7, el nivel de malla con celdas más pequeñas (Scott Wallace, 
2017). Para efectos de este estudio, en la Figura 3, se comparan los resultados de flujo másico reportados 
por Scaggion (2018) y los obtenidos mediante la opción de mallado automático para tres aperturas de 


hadm- 


—9—Scaggion (2018) 
0.04 —E-—Malla 1 


Malla 2 
0,035 —%= Malla 3 
—= Malla 4 
=> 0,03 —98— Malla 5 
99] 
E —+ Malla 6 
20,025 
o 
2 0,02 
Yu 
NS] 
0,015 
pa 
23 0,01 
E 
0,005 
0 
0 1 2 3 4 5 6 7 8 


Apertura de válvula de admisión (mm) 


Figura 3. Comparación de los resultados de flujo másico para hadm= 1,4 y 7 mm obtenidos con los niveles de malla automática 
1,2,3,4,5 y 6 con los reportados por Scaggion (2018). 
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Debido al alto costo computacional que requiere la malla nivel 7 para esta simulación, no fue posible obtener 
los resultados para dicha malla. No obstante, en la gráfica es posible apreciar el comportamiento de las 
variables consideradas para las configuraciones de malla posibles con respecto a la referencia. Se puede 
ver que todos los niveles de malla tienden a converger en hagm= 1 mm a excepción de la malla 1; 
posteriormente se aprecia cierto nivel de divergencia en las magnitudes del flujo másico cuando hadm. = 4 
y 7 mm. Al analizar los tiempos de cómputo para las cinco primeras configuraciones de malla ilustrados en 
la Figura 4, se puede decir que, de todos los niveles de mallas ilustrados, la que más se acerca al 
comportamiento de referencia es la malla 4, cuyo tiempo de cómputo oscila entre los 0,5 y 3,0 minutos. 


10 
9 —9— Malla 1 A Malla 2 
8 Malla 3 —%= Malla 4 


—A= Malla 5 


Tiempo de cómputo (min) 
un 


Apertura de válvula de admisión (mm) 


Figura 4. Tiempos de cómputo para las configuraciones de mallas automáticas 1, 2, 3, 4 y 5. Análisis desarrollado en 
ordenador Lenovo con procesadores Intel Celeron 1.90 Gigahertz y 4 Gigabyte de RAM. 


La configuración de malla seleccionada es la malla 4, la cual cuenta con un promedio de 27,128 celdas 
totales, 12,333 celdas fluidas, 14,795 celdas sólidas y 7,8343 celdas parciales. 


RESULTADOS Y DISCUSIÓN 
Simulación del proceso de admisión 


Con base a los requerimientos generales para una simulación CFD en SWFS, las variables dependientes e 
independientes y las consideraciones antes mencionadas; la Figura 5, expone el diagrama de flujo específico 
propuesto para representar computacionalmente la fluidodinámica al interior del cilindro en estado cuasi 
estacionario durante la carrera de admisión mediante SWFS. 


Ivan Luis Calderón Gutiérrez, Leonardo Redondo Guerra, Robert Macías Naranjo, Fainer Cerpa Olivera, Gail Gutiérrez Ramirez 


6 


Cienciay 
Ingenieria 


Revista Ciencia e Ingeniería, vol. 10, no. 1, eg091885, 2023. 
DOI: https://doi.org/10.5281/zenodo.8091885 


Ly 
Modelo CAD 


Se considera un solo orificio para 
la entrada del fluido al cilindro, 
con: 

Raám = 1,2,3,4,5,6 y 7 mm. 


1. 44 


Fluido de trabajo 
Aire 


Condiciones de pared 
Adiabática con rugosidad (0) 


Condiciones iniciales 
Termodinamicas: 
101325 Pa; 313 K. 
Turbulencia: Por defecto. 
Concentracion: Por defecto. 


WIZARD 


Condiciones de 
frontera 
Condición de 
Presión 
Presión Total 
101325 Pa 


Condición de 
Presión 
Presión Total 
96325 Pa 


Volumen de control 


Entradas 


1. Pared permeable 
Cara/superficie interna 
entrada colector de 
admisión 


Salidas 


1. Pared impermeable 
Caras/superficies del 
Pistón 


Malla 4 automática 
Celdas promedio Totales: 27,128; Fluidas: 
12,333; Solidas: 14,795 y Parciales: 7,843. 


1. Campos de 
velocidad del fluido 


2. Campos de 
presión 


3. Flujo másico en 

pared permeable: 

Entrada Colector 
Admisión 


Figura 5. Diagrama de flujo específico para la simulación CFD del proceso de admisión de un MCI. 


Campos de velocidad y presión en admisión 


En la Figura 6, se muestran los planos de acción que se usarán para ilustrar los campos de velocidad en 
función de fenómenos como Swirl y Tumble, cuyos ejes de rotación son paralelo y perpendicular al cilindro 


respectivamente (El-Adawy et al., 2021). 


(A) Plano Tumble 


(B) Plano Swirl 


Figura 6. Campos de velocidad (A) Tumble y (B) Swirl. Plano Swirl: A, = 15 mm; Tumble B, = 10 mm y B2 = 40 mm. 
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A continuación, la Figura 7, muestra cómo se desarrollan los campos de velocidad desde el plano Tumble 
A, para todo el rango de hadm. Es posible apreciar cómo se intensifican los efectos del fenómeno Tumble 
con la aparición de varios vórtices cuando la apertura de la válvula que permite admisión de aire se vuelve 
mayor, se puede ver como estos remolinos reducen las regiones de velocidad baja ayudando que la 
distribución del aire al interior del cilindro abarque todo el eje longitudinal. 

35.079 7 

31.182 

27.284 

23.386 

19.488 

15.591 

11.693 

7.795 

3.898 

0 


Velocity [m?s] 


35.079 
31.182 4 
27.284 
23.386 
19.488 
15.591 
11.693 
7.795 
3.898 


[1] 
Welocity [m/s] 


(A) Raam = 1 mm (B) Raám = 2 mm (O) haáam = 3 mm (D) haám = 4 mm 


35.079 
Y 31.182 
N 27.284 
23.386 
19.488 
15.591 
11.693 
7,195 
3.898 
0 
MN 


Velocity [mis] 
(PF) haám = 6 mm 


35.079 
31.182 
27.284 
23.386 
19.488 
15.591 
11.693 
7.795 

3.898 


0 
Velocity [m/s] 


(6) Raám =7 mm 


(E) Rhaám = 5 mm 


Figura 7. Distribución del campo de velocidad Tumble en el plano A; para el rango de abertura de la válvula de admisión 
hadm=1,2,3,4,5,6 y 7 mm, 


Los campos de velocidad Swirl en los planos B, y B¿ ubicados a 10 y 40 mm del PMS, se ilustran en la 
Figura 3, teniendo en cuenta la ubicación de la válvula en admisión que se ilustra en higm= 1 mm se hace 
evidente el aumento de la velocidad del fluido de trabajo cerca de las paredes del cilindro en base a la 
aparición de los vórtices. 
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Er [mis] 


35.606 26.609 25.609 
31.650 22.764 22.764 
27.693 19.918 19.918 
23.737 17.073 17.073 
19,781 14227 14227 
15.825 11.382 11.382 
11.869 8.538 8,536 
7.912 5.691 5.591 
3.956 2.845 2,845 
0 o 0 
Velocity [m/s] Velocity [ms] Velocity [rs] 


(4) haam = 1 mm (C) Radm = 3 mm 


35.606 25.509 25.009 
31.650 22.764 22764 
27.693 19.918 19.918 
23.737 17.073 17.073 
19,781 14.227 14.227 
15.826 11.382 11.382 
11.869 8,536 8.536 
7.912 5.691 5.691 
3.956 2.845 2945 
0 0 0 
Velocity [m/s] Velocity [m/s] Velocity [m/s] 


(B) haám =2 mm (D) hadm = 4 mm 


36.606 25.609 
35806 25.609 
31.650 22.764 
27.693 19.918 vic azi 
23737 17.073 paa: als 
19.781 14,227 arar Ha 
AA isa 19781 14227 
11.869 8538 15825 11.382 
7.912 5.691 11.869 8.536 
3.956 2,845 7912 509 
0 0 2.956 2.845 
Velocity [mis] Velocity mis] 0 0 
Velocity [mis] Velocity [ms] 
b. (E) Raám = 5 mm 
h (G) Raam = 7 mm 
35.606 25.609 
31.650 22.764 
27.693 19.918 
2737 17.073 
19.781 14227 
15.825 11.382 
11.869 8.536 
7.912 5.691 
2.956 2.845 
0 0 
Velocity [ms] Velocity ms] 
b, 


(EF) Raam = 6 mm bz 


Velocity [m/s] 


Figura 8. Distribución de los campos de velocidad Swirl en los planos By = 10 mm y B2= 40 mm y para el rango de aberturas 
de la válvula de admisión hadm= 1,2,3,4,5, 6 y 7 mm. 


Los fenómenos Tumble y Swirl antes ilustrados, se originan durante la carrera de admisión y sus efectos se 
hacen evidentes al interior del cilindro influyendo en aspectos como formación de la mezcla, guiado de la 
mezcla y en las condiciones de turbulencia del fluido de trabajo; es importante resaltar que los efectos Swirl 
tienden a incrementar en las proximidades del Punto Muerto Inferior (PMI) (Vázquez áú De la Morena, 
2019), por lo cual en el campo B, se tienen unas características de flujo más ordenadas. La Figura 9, 
permite apreciar los campos de presión para todo el rango hagm con base a la diferencia de presiones que 
se configuró como condición de frontera. 
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101330.50 
100694.72 
100058.94 
9942317 
98787.39 
98151.61 
97515.83 
96880.06 
96244 28 


95608.50 
Pressure [Pal (4) Radám = 1 mm 


(C) haám = 3 mm (D) haam = 4 mm 


101330.50 
100694.72 
100058.94 
9942317 
98787.39 
98151.61 
97515.83 
96880.06 
9624428 


95608.50 
Pressure [Pa] 


(E) haám = 5 mm (E) haám = 6 mm (G) haám = 7 mm 


Figura 9. Comportamiento de los campos de presión en el rango de aberturas de válvula de hadm = 1, 2,3, 4,5, 6 y 7 mm. 


Resultados cuantitativos en admisión 


Considerando que el flujo másico admitido se seleccionó como una variable dependiente en la simulación 
de la carrera de admisión, es decir como un objetivo de la simulación, en la Tabla 2 se muestran las 
magnitudes calculadas para todo el rango de aperturas hadm. 


Tabla 2. Magnitudes del flujo másico admitido para todo el rango de aperturaS hadm- 


Rango aperturas hadm | Colector admisión: flujo másico 
(mm) (kg.s”) 

0,0072 

0,0131 

0,0194 

0,0248 

0,0282 

0,0314 

0,0339 


N|OMO|U|Aa|w| nj] — 
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Se evidencia un incremento del flujo másico a medida que hadm se hace más grande, este aumento se da 
porque el área mínima de paso del flujo de aire crece cuando la apertura de válvula hadm aumenta y la 
limitación geométrica que limita el paso del flujo de aire disminuye, permitiendo de esta manera dicho 
incremento. 


Simulación del proceso de escape 


Para la simulación de la carrera de escape se ajustó el flujo de trabajo ¡ilustrado en la Figura 10, en dicho 
flujo de trabajo se pueden apreciar las condiciones iniciales, de frontera, consideraciones y objetivos de 
simulación anteriormente mencionados para este escenario en la subsección condiciones iniciales y de 
frontera. 


Modelo CAD 
Se considera un solo orificio para 


colector admisión con (Raqm= a Volumen de control 


Condiciones de 


4 mm) y colector escape con 


(Resc= 1,4 y 10 mm) frontera Entradas 
Condición de Flujo 1. Pared permeable 1. Campos de velocidad 
my Entrada de Flujo Cara/superficie interna : 
0,0248 Kg/s entrada colector admisión dun), 2. Campos de presión 


Fluido de trabajo 
Aire 


Condiciones de pared m) 
] 
] 
] 
] 


Salidas 


1. Pared permeable 
Cara/superficie interna 
salida colector escape 


| 
l 3. Velocidad en pared 
I permeable: Entrada 
' colector admisión y 

' salida colector escape. 


Condición de Presión 
Presión Total 
101325 Pa 


Goals 


Adiabática con rugosidad (0) 


WIZARD 


4. Flujo másico en 
pared permeable: 


Condiciones iniciales 
Termodinamicas: 
101325 Pa; 313 K. 
Turbulencia: Por defecto. 
Concentracion: Por defecto. 


12,333; Solidas: 14,795 y Parciales: 7,843. 


| 
JS | 
Malla 4 automática l Salida colector escape. 
Celdas promedio Totales: 27,128; Fluidas: I 
| 
| 


Figura 10. Diagrama de flujo específico para la simulación CFD del proceso de escape de un MCI, 


Campos de velocidad y presión en escape 


Los campos de velocidad se ilustran en la Figura 11, estos muestran la admisión y posterior escape del 
flujo de aire mientras al mantener la abertura en admisión hadm = 4 mm y variar las aberturas de escape 
según el rango hes.= 1,4 y 10 mm. Con base a estas condiciones es posible observar las características 
del flujo de aire mientras transita por el modelo CAD y sus respectivas variaciones. La Figura 12, muestra 
los campos de presión, con base en estos es posible identificar el proceso de transferencia de masa de un 
punto a otro considerando que dicho proceso se presenta de zonas con altas presiones a las zonas de 
bajas presiones hasta llegar a un equilibrio entre las presiones en cuestión (Cengel et al., 2019). 
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90.036 
80.032 
70.028 
50.024 
50.020 
40.016 
30.012 
20.008 
10.004 
0 

Velocity [mis] 


(A) hos. = 1 mm (B) hos. = ¿mm 


(C) hos. = 10 mm 


Figura 11. Campos de velocidad del flujo de aire para las condiciones geométricas hadm = 4 mm y hesc = 1,4 y 10 mm. 


117792.96 
115499 60 
113206.25 
11091289 
108619.54 
106326.18 
10403282 
101739.47 
99446.11 


97152.75 
Pressure [Pal 


(A) Raso = 1 mm (B) hos. = 4 mm 


109843.03 
10860313 
107363.24 
106123.34 
10498344 
103643.54 
102403.64 
101163.74 
9992384 
93683.95 
Pressure [Pa] 


(C) Res. = 10 mm 


Figura 12. Campos de presión para las condiciones geométricas hadm = 4 MM y hese = 1,4 y 10 mm. 
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Resultados cuantitativos en escape 


En la Figura 13. se ilustra el comportamiento de las velocidades de admisión y escape con base a las 
magnitudes calculadas de esta variable en función del rango de apertura hesc. 


20 85 
E Ta 84 
E” E 
3 16 2 83 
«L 3 82 
E 14 1 Se 
[e] 
a 12 3 81 
a E 
ES 
y 
Z 10 E 80 
3 2 
2 8 vo 79 
E > 
6 78 
0.112.344 56 7 8 9 10 11 0.112.3.4.556 7 8 9 10 11 
Aberturas válvula de escape (mm) Aberturas válvula de escape (mm) 


Figura 13. Velocidades de admisión y escape en función de las aberturas de escape hese = 1,4 y 10 mm. 


La velocidad de admisión incrementa a medida que hesc aumenta, se puede decir que a partir de has. = 4 
mm, la velocidad de admisión muestra un comportamiento constante, lo cual se debe a que el área de paso 
de la geometría en la zona de admisión con la apertura de válvula hagm = 4 mm no permite admitir más 
alre representado de esta manera una limitante geométrica. En cuanto a la velocidad de escape esta muestra 
el comportamiento contrario y se debe a que el área de paso del flujo de aire en la zona de escape aumenta, 
representado una caída de velocidad a medida que has. disminuye. 


Finalmente, en la Tabla 3. se muestran las magnitudes del flujo másico calculadas en la salida del colector 
de escape. 


Tabla 3. Magnitudes del flujo másico medido a la salida del colector de escape en función del rango de apertura hesc = 1,4 y 
10 mm. 


Rango apertura hosc Colector escape: flujo 
(mm) másico (kg.s”*) 
1 -0,0248 
4 -0,0248 
10 -0,0248 
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Validación de resultados 


Se graficaron los resultados del flujo másico reportados por la literatura referencia y los obtenidos en este 
estudio para el rango de aberturas de admisión hadm = 1, 2,3, 4,5, 6 y 7 mm lo cual se ilustra en la Figura 
14, 


—$—Scaggion (2018) 
—lERSWFS 


Flujo másico (kg.s"!) 


0 l 2 3 + 5 6 7 8 
Aperturas válvula de admisión (mm) 


Figura 14. Comparación de los resultados obtenidos mediante SWFS y los reportados por Scaggion (2018). 
Para las aperturas hadm = 1 y 2 mm se puede ver que los resultados convergen en cierto grado y a partir 
de hadm = 3 mm divergen aumentando dicho nivel hasta llegar a hadm = 7 mm. Sin embargo, esta 


discrepancia se debe a las diferencias geométricas existentes entre el modelo CAD usado por Scaggion 
(2018) y el usado en este estudio, las cuales se listan en la Tabla 4. 


Tabla 4. Diferencias geométricas entre los modelos CAD usados. 


Parámetros Scaggion (2018) SWFS 
Ángulo de inclinación de válvula admisión 11,509 19,90% 
Ángulo de inclinación válvula escape 12,009 19,909 
Radio válvula admisión 13,75 mm 13,65 mm 
Radio válvula escape 11,00 mm 11,6 mm 
Abertura máx. válvula admisión 8,3 mm 7,0 mm 
Abertura máx. válvula escape 7,2 mm 10,0 mm 
Radio del cilindro 33,50 mm 35,95 


Para la validación del proceso de escape se aplicó el principio de conservación de masa "Min — Mout = 
dmyc/dt" reportado por (Cengel et al., 2019), concluyendo que el flujo másico que entra al cilindro fue el 
mismo que salió para todo el rango hes. considerado; por lo tanto, se puede decir que esta geometría no 
retiene gases de escape. 
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CONCLUSIONES 


+ 


Fue posible representar la fluidodinámica del aire al interior del cilindro en términos de campos de velocidad 
y presión en condiciones cuasi estacionarias para todo el rango considerado de aperturas hadm: 
especialmente la visualización de fenómenos Tumble y Swirl los cuales pueden favorecer significativamente 
a la homogeneización del fluido de trabajo en todo el cilindro durante la carrera de admisión e igualmente 
al proceso de transferencia de masa. 


La fluidodinámica presente en la carrera de escape indica que el modelo CAD usado en este estudio es 
viable para análisis futuros debido a que no retiene gases de escape en todo el rango de aperturas hesc. 


La validación de resultados pone en evidencia que las metodologías permiten la representación de la 
fluidodinámica del flujo de aire en el cilindro con un alto nivel de convergencia. Lo cual indica que SWFS es 
una herramienta CFD con capacidad para el estudio de factores claves como las cualidades de la formación 
de la mezcla aire/combustible y el movimiento de esta mezcla al interior del cilindro en el funcionamiento 
de los motores y que podría usarse en el estudio de procesos innovadores como el Blending con la inclusión 
de combustibles renovables en los MCI. 
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